Quantitative trait loci and genomic prediction for grain sugar and mineral concentrations of cowpea [Vigna unguiculata (L.) Walp.].

Development of high yielding cowpea varieties coupled with good taste and rich in essential minerals can promote consumption and thus nutrition and profitability. The sweet taste of cowpea grain is determined by its sugar content, which comprises mainly sucrose and galacto-oligosaccharides (GOS) including raffinose and stachyose. However, GOS are indigestible and their fermentation in the colon can produce excess intestinal gas, causing undesirable bloating and flatulence. In this study, we aimed to examine variation in grain sugar and mineral concentrations, then map quantitative trait loci (QTLs) and estimate genomic-prediction (GP) accuracies for possible application in breeding. Grain samples were collected from a multi-parent advanced generation intercross (MAGIC) population grown in California during 2016-2017. Grain sugars were assayed using high-performance liquid chromatography. Grain minerals were determined by inductively coupled plasma-optical emission spectrometry and combustion. Considerable variation was observed for sucrose (0.6-6.9%) and stachyose (2.3-8.4%). Major QTLs for sucrose (QSuc.vu-1.1), stachyose (QSta.vu-7.1), copper (QCu.vu-1.1) and manganese (QMn.vu-5.1) were identified. Allelic effects of major sugar QTLs were validated using the MAGIC grain samples grown in West Africa in 2017. GP accuracies for minerals were moderate (0.4-0.58). These findings help guide future breeding efforts to develop mineral-rich cowpea varieties with desirable sugar content.

A novel natural variation in the promoter of GmCHX1 regulates conditional gene expression to improve salt tolerance in soybean.

Identification and characterization of soybean germplasm and gene(s)/allele(s) for salt tolerance is an effective way to develop improved varieties for saline soils. Previous studies identified GmCHX1 (Glyma03g32900) as a major salt tolerance gene in soybean, and two main functional variations were found in the promoter region (148/150 bp insertion) and the third exon with a retrotransposon insertion (3.78 kb). In the current study, we identified four salt-tolerant soybean lines, including PI 483460B (Glycine soja), carrying the previously identified salt-sensitive variations at GmCHX1, suggesting new gene(s) or new functional allele(s) of GmCHX1 in these soybean lines. Subsequently, we conducted quantitative trait locus (QTL) mapping in a recombinant-inbred line population (Williams 82 (salt-sensitive) × PI 483460B) to identify the new salt tolerance loci/alleles. A new locus, qSalt_Gm18, was mapped on chromosome 18 associated with leaf scorch score. Another major QTL, qSalt_Gm03, was identified to be associated with chlorophyll content ratio and leaf scorch score in the same chromosomal region of GmCHX1 on chromosome 3. Novel variations in a STRE (stress response element) cis-element in the promoter region of GmCHX1 were found to regulate the salt-inducible expression of the gene in these four newly identified salt-tolerant lines including PI 483460B. This new allele of GmCHX1 with salt-inducible expression pattern provides an energy cost efficient (conditional gene expression) strategy to protect soybean yield in saline soils without yield penalty under non-stress conditions. Our results suggest that there might be no other major salt tolerance locus similar to GmCHX1 in soybean germplasm, and further improvement of salt tolerance in soybean may rely on gene-editing techniques instead of looking for natural variations.

Novel genetic resources associated with sucrose and stachyose content through genome-wide association study in soybean (Glycine max (L.) Merr.).

The nutritional value of soybean [Glycine max (L.) Merr.] for animals is influenced by soluble carbohydrates, such as sucrose and stachyose. Although sucrose is nutritionally desirable, stachyose is an antinutrient causing diarrhea and flatulence in non-ruminant animals. We conducted a genome-wide association study of 220 soybean accessions using 21,317 single nucleotide polymorphisms (SNPs) from the SoySNP50K iSelect Beadchip data to identify significant SNPs associated with sucrose and stachyose content. Seven significant SNPs were identified for sucrose content across chromosomes (Chrs.) 2, 8, 12, 17, and 20, while thirteen significant SNPs were identified for stachyose content across Chrs. 2, 5, 8, 9, 10, 13, 14, and 15. Among those significant SNPs, three sucrose-related SNPs on Chrs. 8 and 17 were novel, while twelve stachyose-related SNPs on Chrs. 2, 5, 8, 9, 10, 13, 14, and 15 were novel. Based on Phytozome, STRING, and GO annotation, 17 and 24 candidate genes for sucrose and stachyose content, respectively, were highly associated with the carbohydrate metabolic pathway. Among these, the publicly available RNA-seq Atlas database highlighted four candidate genes associated with sucrose (Glyma.08g361200 and Glyma.17g258100) and stachyose (Glyma.05g025300 and Glyma.13g077900) content, which had higher gene expression levels in developing seed and multiple parts of the soybean plant. The results of this study will extend knowledge of the molecular mechanism and genetic basis underlying sucrose and stachyose content in soybean seed. Furthermore, the novel candidate genes and SNPs can be valuable genetic resources that soybean breeders may utilize to modify carbohydrate profiles for animal and human usage.

Genomic analysis and characterization of new loci associated with seed protein and oil content in soybeans.

Breeding for increased protein without a reduction in oil content in soybeans [Glycine max (L.) Merr.] is a challenge for soybean breeders but an expected goal. Many efforts have been made to develop new soybean varieties with high yield in combination with desirable protein and/or oil traits. An elite line, R05-1415, was reported to be high yielding, high protein, and low oil. Several significant quantitative trait loci (QTL) for protein and oil were reported in this line, but many of them were unstable across environments or genetic backgrounds. Thus, a new study under multiple field environments using the Infinium BARCSoySNP6K BeadChips was conducted to detect and confirm stable genomic loci for these traits. Genetic analyses consistently detected a single major genomic locus conveying these two traits with remarkably high phenotypic variation explained (R2 ), varying between 24.2% and 43.5%. This new genomic locus is located between 25.0 and 26.7 Mb, distant from the previously reported QTL and did not overlap with other commonly reported QTL and the recently cloned gene Glyma.20G085100. Homolog analysis indicated that this QTL did not result from the paracentric chromosome inversion with an adjacent genomic fragment that harbors the reported QTL. The pleiotropic effect of this QTL could be a challenge for improving protein and oil simultaneously; however, a further study of four candidate genes with significant expressions in the seed developmental stages coupled with haplotype analysis may be able to pinpoint causative genes. The functionality and roles of these genes can be determined and characterized, which lay a solid foundation for the improvement of protein and oil content in soybeans.

Near-gapless genome assemblies of Williams 82 and Lee cultivars for accelerating global soybean research.

Complete, gapless telomere-to-telomere chromosome assemblies are a prerequisite for comprehensively investigating the architecture of complex regions, like centromeres or telomeres and removing uncertainties in the order, spacing, and orientation of genes. Using complementary genomics technologies and assembly algorithms, we developed highly contiguous, nearly gapless, genome assemblies for two economically important soybean [Glycine max (L.) Merr] cultivars (Williams 82 and Lee). The centromeres were distinctly annotated on all the chromosomes of both assemblies. We further found that the canonical telomeric repeats were present at the telomeres of all chromosomes of both Williams 82 and Lee genomes. A total of 10 chromosomes in Williams 82 and eight in Lee were entirely reconstructed in single contigs without any gap. Using the combination of ab initio prediction, protein homology, and transcriptome evidence, we identified 58,287 and 56,725 protein-coding genes in Williams 82 and Lee, respectively. The genome assemblies and annotations will serve as a valuable resource for studying soybean genomics and genetics and accelerating soybean improvement.

Classification Methods and Identification of Reniform Nematode Resistance in Known Soybean Cyst Nematode-Resistant Soybean Genotypes.

Plant parasitic nematodes are a major yield-limiting factor of soybean in the United States and Canada. It has been indicated that soybean cyst nematode (SCN; Heterodera glycines Ichinohe) and reniform nematode (RN; Rotylenchulus reniformis Linford and Oliveira) resistance could be genetically related. For many years, fragmentary data have shown this relationship. This report evaluates RN reproduction on 418 plant introductions (PIs) selected from the U.S. Department of Agriculture Soybean Germplasm Collection with reported SCN resistance. The germplasm was divided into two tests of 214 PIs reported as resistant and 204 PIs reported as moderately resistant to SCN. The defining and reporting of RN resistance changed several times in the last 30 years, causing inconsistencies in RN resistance classification among multiple experiments. Comparison of four RN resistance classification methods was performed: (i) ≤10% as compared with the susceptible check, (ii) using normalized reproduction index (RI) values, and using (iii) transformed data log10(x), and (iv) transformed data log10(x + 1) in an optimal univariate k-means clustering analysis. The method of transformed data log10(x) was selected as the most accurate for classification of RN resistance. Among 418 PIs with reported SCN resistance, the log10(x) method grouped 59 PIs (15%) as resistant and 130 PIs (31%) as moderately resistant to RN. Genotyping of a subset of the most resistant PIs to both nematode species revealed their strong correlation with rhg1-a allele. This research identified genotypes with resistance to two nematode species and potential new sources of RN resistance that could be valuable to breeders in developing resistant cultivars.

Sequencing the USDA core soybean collection reveals gene loss during domestication and breeding.

The gene content of plants varies between individuals of the same species due to gene presence/absence variation, and selection can alter the frequency of specific genes in a population. Selection during domestication and breeding will modify the genomic landscape, though the nature of these modifications is only understood for specific genes or on a more general level (e.g., by a loss of genetic diversity). Here we have assembled and analyzed a soybean (Glycine spp.) pangenome representing more than 1,000 soybean accessions derived from the USDA Soybean Germplasm Collection, including both wild and cultivated lineages, to assess genomewide changes in gene and allele frequency during domestication and breeding. We identified 3,765 genes that are absent from the Lee reference genome assembly and assessed the presence/absence of all genes across this population. In addition to a loss of genetic diversity, we found a significant reduction in the average number of protein-coding genes per individual during domestication and subsequent breeding, though with some genes and allelic variants increasing in frequency associated with selection for agronomic traits. This analysis provides a genomic perspective of domestication and breeding in this important oilseed crop.

Mapping of partial resistance to Phytophthora sojae in soybean PIs using whole-genome sequencing reveals a major QTL.

In the last decade, more than 70 quantitative trait loci (QTL) related to soybean [Glycine max (L.) Merr.] partial resistance (PR) against Phytophthora sojae have been identified by genome-wide association studies (GWAS). However, most of them have either a minor effect on the resistance level or are specific to a single phenotypic variable or one isolate, thereby limiting their use in breeding programs. In this study, we have used an analytical approach combining (a) the phenotypic characterization of a diverse panel of 357 soybean accessions for resistance to P. sojae captured through a single variable, corrected dry weight; (b) a new hydroponic assay allowing the inoculation of a combination of P. sojae isolates covering the spectrum of commercially relevant Rps genes; and (c) exhaustive genotyping through whole-genome resequencing (WGS). This led to the identification of a novel P. sojae resistance QTL with a relatively major effect compared with the previously reported QTL. The QTL interval, spanning ∼500 kb on chromosome (Chr) 15, does not colocalize with previously reported QTL for P. sojae resistance. Plants carrying the favorable allele at this QTL were 60% more resistant. Eight genes were found to reside in the linkage disequilibrium (LD) block containing the peak single-nucleotide polymorphism (SNP) including Glyma.15G217100, which encodes a major latex protein (MLP)-like protein, with a functional annotation related to pathogen resistance. Expression analysis of Glyma.15G217100 indicated that it was nearly eight times more highly expressed in a group of plant introductions (PIs) carrying the resistant (R) allele compared with those carrying the susceptible (S) allele within a short period after inoculation. These results offer new and valuable options to develop improved soybean cultivars with broad resistance to P. sojae through marker-assisted selection.

Genome-wide identification and analysis of soybean acyl-ACP thioesterase gene family reveals the role of GmFAT to improve fatty acid composition in soybean seed.

KEY MESSAGE: Soybean acyl-ACP thioesterase gene family have been characterized; GmFATA1A mutants were discovered to confer high oleic acid, while GmFATB mutants presented low palmitic and high oleic acid seed content. Soybean oil stability and quality are primarily determined by the relative proportions of saturated versus unsaturated fatty acids. Commodity soybean typically contains 11% palmitic acid, as the primary saturated fatty acids. Reducing palmitic acid content is the principal approach to minimize the levels of saturated fatty acids in soybean. Though high palmitic acid enhances oxidative stability of soybean oil, it is negatively correlated with oil and oleic acid content and can cause coronary heart diseases for humans. For plants, acyl-acyl carrier protein (ACP) thioesterases (TEs) are a group of enzymes to hydrolyze acyl group and release free fatty acid from plastid. Among them, GmFATB1A has become the main target to genetically reduce the palmitic acid content in soybean. However, the role of members in soybean acyl-ACP thioesterase gene family is largely unknown. In this study, we characterized two classes of TEs, GmFATA, and GmFATB in soybean. We also denominated two GmFATA members and discovered six additional members that belong to GmFATB gene family through phylogenetic, syntenic, and in silico analysis. Using TILLING-by-Sequencing+, we identified an allelic series of mutations in five soybean acyl-ACP thioesterase genes, including GmFATA1A, GmFATB1A, GmFATB1B, GmFATB2A, and GmFATB2B. Additionally, we discovered mutations at GmFATA1A to confer high oleic acid (up to 34.5%) content, while mutations at GmFATB presented low palmitic acid (as low as 5.6%) and high oleic acid (up to 36.5%) phenotypes. The obtained soybean mutants with altered fatty acid content can be used in soybean breeding program for improving soybean oil composition traits.

TILLING-by-Sequencing+ Reveals the Role of Novel Fatty Acid Desaturases (GmFAD2-2s) in Increasing Soybean Seed Oleic Acid Content.

Soybean is the second largest source of oil worldwide. Developing soybean varieties with high levels of oleic acid is a primary goal of the soybean breeders and industry. Edible oils containing high level of oleic acid and low level of linoleic acid are considered with higher oxidative stability and can be used as a natural antioxidant in food stability. All developed high oleic acid soybeans carry two alleles; GmFAD2-1A and GmFAD2-1B. However, when planted in cold soil, a possible reduction in seed germination was reported when high seed oleic acid derived from GmFAD2-1 alleles were used. Besides the soybean fatty acid desaturase (GmFAD2-1) subfamily, the GmFAD2-2 subfamily is composed of five members, including GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E. Segmental duplication of GmFAD2-1A/GmFAD2-1B, GmFAD2-2A/GmFAD2-2C, GmFAD2-2A/GmFAD2-2D, and GmFAD2-2D/GmFAD2-2C have occurred about 10.65, 27.04, 100.81, and 106.55 Mya, respectively. Using TILLING-by-Sequencing+ technology, we successfully identified 12, 8, 10, 9, and 19 EMS mutants at the GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E genes, respectively. Functional analyses of newly identified mutants revealed unprecedented role of the five GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E members in controlling the seed oleic acid content. Most importantly, unlike GmFAD2-1 members, subcellular localization revealed that members of the GmFAD2-2 subfamily showed a cytoplasmic localization, which may suggest the presence of an alternative fatty acid desaturase pathway in soybean for converting oleic acid content without substantially altering the traditional plastidial/ER fatty acid production.

Identification and characterization of novel QTL conferring internal detoxification of aluminium in soybean.

Aluminium (Al) toxicity inhibits soybean root growth, leading to insufficient water and nutrient uptake. Two soybean lines ('Magellan' and PI 567731) were identified differing in Al tolerance, as determined by primary root length ratio, total root length ratio, and root tip number ratio under Al stress. Serious root necrosis was observed in PI 567731, but not in Magellan under Al stress. An F8 recombinant inbred line population derived from a cross between Magellan and PI 567731 was used to map the quantitative trait loci (QTL) for Al tolerance. Three QTL on chromosomes 3, 13, and 20, with tolerant alleles from Magellan, were identified. qAl_Gm13 and qAl_Gm20 explained large phenotypic variations (13-27%) and helped maintain root elongation and initiation under Al stress. In addition, qAl_Gm13 and qAl_Gm20 were confirmed in near-isogenic backgrounds and were identified to epistatically regulate Al tolerance via internal detoxification instead of Al3+ exclusion. Phylogenetic and pedigree analysis identified the tolerant alleles of both loci derived from the US ancestral line, A.K.[FC30761], originally from China. Our results provide novel genetic resources for breeding Al-tolerant soybean and suggest that internal detoxification contributes to soybean tolerance to excessive soil Al.

Dissecting nematode resistance regions in soybean revealed pleiotropic effect of soybean cyst and reniform nematode resistance genes.

Reniform nematode (RN, Rotylenchulus reniformis Linford & Oliveira) has emerged as one of the most important plant parasitic nematodes of soybean [Glycine max (L.) Merr.]. Planting resistant varieties is the most effective strategy for nematode management. The objective of this study was to identify quantitative trait loci (QTL) for RN resistance in an exotic soybean line, PI 438489B, using two linkage maps constructed from the Universal Soybean Linkage Panel (USLP 1.0) and next-generation whole-genome resequencing (WGRS) technology. Two QTL controlling RN resistance were identified-the soybean cyst nematode (SCN, Heterodera glycines) resistance gene GmSNAP18 at the rhg1 locus and its paralog GmSNAP11. Strong association between resistant phenotype and haplotypes of the GmSNAP11 and GmSNAP18 was observed. The results indicated that GmSNAP11 possibly could have epistatic effect on GmSNAP18, or vice versa, with the presence of a significant correlation in RN resistance of rhg1-a GmSNAP18 vs. rhg1-b GmSNAP18. Most importantly, our preliminary data suggested that GmSNAP18 and GmSNAP11 proteins physically interact in planta, suggesting that they belong to the same pathway for resistance. Unlike GmSNAP18, no indication of GmSNAP11 copy number variation was found. Moreover, gene-based single nucleotide polymorphism (SNP) markers were developed for rapid detection of RN or SCN resistance at these loci. Our analysis substantiates synergic interaction between GmSNAP11 and GmSNAP18 genes and confirms their roles in RN as well as SCN resistance. These results could contribute to a better understanding of evolution and subfunctionalization of genes conferring resistance to multiple nematode species and provide a framework for further investigations.

Genetic characterization of qSCN10 from an exotic soybean accession PI 567516C reveals a novel source conferring broad-spectrum resistance to soybean cyst nematode.

KEY MESSAGE: The qSCN10 locus with broad-spectrum SCN resistance was fine-mapped to a 379-kb region on chromosome 10 in soybean accession PI 567516C. Candidate genes and potential application benefits of this locus were discussed. Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is one of the most devastating pests of soybean, causing significant yield losses worldwide every year. Genetic resistance has been the major strategy to control this pest. However, the overuse of the same genetic resistance derived primarily from PI 88788 has led to the genetic shifts in nematode populations and resulted in the reduced effectiveness in soybean resistance to SCN. Therefore, novel genetic resistance resources, especially those with broad-spectrum resistance, are needed to develop new resistant cultivars to cope with the genetic shifts in nematode populations. In this study, a quantitative trait locus (QTL) qSCN10 previously identified from a soybean landrace PI 567516C was confirmed to confer resistance to multiple SCN HG Types. This QTL was further fine-mapped to a 379-kb region. There are 51 genes in this region. Four of them are defense-related and were regulated by SCN infection, suggesting their potential role in mediating resistance to SCN. The phylogenetic and haplotype analyses of qSCN10 as well as other information indicate that this locus is different from other reported resistance QTL or genes. There was no yield drag or other unfavorable traits associated with this QTL when near-isogenic lines with and without qSCN10 were tested in a SCN-free field. Therefore, our study not only provides further insight into the genetic basis of soybean resistance to SCN, but also identifies a novel genetic resistance resource for breeding soybean for durable, broad-spectrum resistance to this pest.

Fine-mapping and characterization of qSCN18, a novel QTL controlling soybean cyst nematode resistance in PI 567516C.

KEY MESSAGE: The qSCN18 QTL from PI 56756C was confirmed and fine-mapped to improve soybean resistance to the SCN population HG Type 2.5.7 using near-isogenic lines carrying recombination crossovers within the QTL region. The QTL underlying resistance was fine-mapped to a 166-Kbp region on chromosome 18, and the candidate genes were selected based on genomic analyses. Soybean cyst nematode (SCN, Heterodera glycines, Ichinohe) is the most devastating pathogen of soybean. Understanding the genetic basis of SCN resistance is crucial for managing this parasite in the field. Two major loci, rhg1 and Rhg4, were previously characterized as valuable resources for SCN resistance. However, their continuous use has caused shifts in the virulence of SCN populations, which can overcome the resistance conferred by these two major loci. Reduced effectiveness became a major concern in the soybean industry due to continuous use of rhg1 for decades. Thus, it is imperative to identify sources of SCN resistance for durable SCN management. A novel QTL qSCN18 was identified in PI567516C. To fine-map qSCN18 and identify resistance genes, a large backcross population was developed. Nineteen near-isogenic lines (NILs) carrying recombination crossovers within the QTL region were identified. The first phase of fine-mapping narrowed the QTL region to 549-Kbp, whereas the second phase confined the region to 166-Kbp containing 23 genes. Two flanking markers, MK-1 and MK-6, were developed and validated to detect the presence of the qSCN18 resistance allele. Haplotype analysis clustered the fine-mapped qSCN18 region from PI 567516C with the cqSCN-007 locus previously mapped in the wild soybean accession PI 468916. The NILs were developed to further characterize the causal gene(s) harbored in this QTL. This study also confirmed the previously identified qSCN18. The results will facilitate marker-assisted selection (MAS) introducing the qSCN18 locus from PI 567516C into high-yielding soybean cultivars with durable resistance to SCN.

Mapping Quantitative Trait Loci for Soybean Seedling Shoot and Root Architecture Traits in an Inter-Specific Genetic Population.

Wild soybean species (Glycine soja Siebold & Zucc.) comprise a unique resource to widen the genetic base of cultivated soybean [Glycine max (L.) Merr.] for various agronomic traits. An inter-specific mapping population derived from a cross of cultivar Williams 82 and PI 483460B, a wild soybean accession, was utilized for genetic characterization of root architecture traits. The objectives of this study were to identify and characterize quantitative trait loci (QTL) for seedling shoot and root architecture traits, as well as to determine additive/epistatic interaction effects of identified QTLs. A total of 16,469 single nucleotide polymorphisms (SNPs) developed for the Illumina beadchip genotyping platform were used to construct a high resolution genetic linkage map. Among the 11 putative QTLs identified, two significant QTLs on chromosome 7 were determined to be associated with total root length (RL) and root surface area (RSA) with favorable alleles from the wild soybean parent. These seedling root traits, RL (BARC_020495_04641 ~ BARC_023101_03769) and RSA (SNP02285 ~ SNP18129_Magellan), could be potential targets for introgression into cultivated soybean background to improve both tap and lateral roots. The RL QTL region harbors four candidate genes with higher expression in root tissues: Phosphofructokinase (Glyma.07g126400), Snf7 protein (Glyma.07g127300), unknown functional gene (Glyma.07g127900), and Leucine Rich-Repeat protein (Glyma.07g127100). The novel alleles inherited from the wild soybean accession could be used as molecular markers to improve root system architecture and productivity in elite soybean lines.

Whole-genome resequencing identifies quantitative trait loci associated with mycorrhizal colonization of soybean.

KEY MESSAGE: A whole-genome resequencing-derived SNP dataset identified six quantitative trait loci (QTL) significantly associated with colonization of soybean by an arbuscular mycorrhizal fungus (Rhizophagus intraradices). Candidate genes identified in these QTL regions include homologs to known nodulin protein families and other symbiosis-specific genes. Arbuscular mycorrhizal fungi (AMF) form associations with over 80% of all terrestrial plant species and assist their host plants by increasing their nutrient uptake, drought tolerance, and resilience against pathogens and pests. Genotypic variation of crop plants to AMF colonization has been identified in crops, including soybean; however, the genetics controlling levels of AMF colonization in soybean are unknown. The overall goal of our study was to identify genomic regions associated with mycorrhizal colonization in soybean using genome-wide association analysis. A diverse panel of 350 exotic soybean genotypes inoculated with Rhizophagus intraradices were microscopically evaluated for root colonization using a modified gridline intersect method. Root colonization differed significantly (P < 0.001) among genotypes and ranged from 11 to 70%. A whole-genome resequencing-derived SNP dataset identified six quantitative trait loci (QTL) significantly associated with R. intraradices colonization that explained 24% of the phenotypic variance. Candidate genes identified in these QTL regions include homologs to known nodulin protein families and other symbiosis-specific genes. The results showed there was a significant genetic component to the level of colonization by R. intraradices in soybean. This information may be useful in the development of AMF-sensitive soybean cultivars to enhance nutrient uptake, drought tolerance, and disease resistance in the crop.

Assessment of Phenotypic Variations and Correlation among Seed Composition Traits in Mutagenized Soybean Populations.

Soybean [Glycine max (L.) Merr.] seed is a valuable source of protein and oil worldwide. Traditionally, the natural variations were heavily used in conventional soybean breeding programs to select desired traits. However, traditional plant breeding is encumbered with low frequencies of spontaneous mutations. In mutation breeding, genetic variations from induced mutations provide abundant sources of alterations in important soybean traits; this facilitated the development of soybean germplasm with modified seed composition traits to meet the different needs of end users. In this study, a total of 2366 'Forrest'-derived M2 families were developed for both forward and reverse genetic studies. A subset of 881 M3 families was forward genetically screened to measure the contents of protein, oil, carbohydrates, and fatty acids. A total of 14 mutants were identified to have stable seed composition phenotypes observed in both M3 and M4 generations. Correlation analyses have been conducted among ten seed composition traits and compared to a collection of 103 soybean germplasms. Mainly, ethyl methanesulfonate (EMS) mutagenesis had a strong impact on the seed-composition correlation that was observed among the 103 soybean germplasms, which offers multiple benefits for the soybean farmers and industry to breed for desired multiple seed phenotypes.

Identification of new loci for salt tolerance in soybean by high-resolution genome-wide association mapping.

BACKGROUND: Salinity is an abiotic stress that negatively affects soybean [Glycine max (L.) Merr.] seed yield. Although a major gene for salt tolerance was identified and consistently mapped to chromosome (Chr.) 3 by linkage mapping studies, it does not fully explain genetic variability for tolerance in soybean germplasm. In this study, a genome-wide association study (GWAS) was performed to map genomic regions for salt tolerance in a diverse panel of 305 soybean accessions using a single nucleotide polymorphism (SNP) dataset derived from the SoySNP50K iSelect BeadChip. A second GWAS was also conducted in a subset of 234 accessions using another 3.7 M SNP dataset derived from a whole-genome resequencing (WGRS) study. In addition, three gene-based markers (GBM) of the known gene, Glyma03g32900, on Chr. 3 were also integrated into the two datasets. Salt tolerance among soybean lines was evaluated by leaf scorch score (LSS), chlorophyll content ratio (CCR), leaf sodium content (LSC), and leaf chloride content (LCC). RESULTS: For both association studies, a major locus for salt tolerance on Chr. 3 was confirmed by a number of significant SNPs, of which three gene-based SNP markers, Salt-20, Salt14056 and Salt11655, had the highest association with all four traits studied. Also, additional genomic regions on Chrs. 1, 8, and 18 were found to be associated with various traits measured in the second GWAS using the WGRS-derived SNP dataset. CONCLUSIONS: A region identified on Chr. 8 was identified to be associated with all four traits and predicted as a new minor locus for salt tolerance in soybean. The candidate genes harbored in this minor locus may help reveal the molecular mechanism involved in salt tolerance and to improve tolerance in soybean cultivars. The significant SNPs will be useful for marker-assisted selection for salt tolerance in soybean breeding programs.

Whole-genome re-sequencing reveals the impact of the interaction of copy number variants of the rhg1 and Rhg4 genes on broad-based resistance to soybean cyst nematode.

Soybean cyst nematode (SCN) is the most devastating plant-parasitic nematode. Most commercial soybean varieties with SCN resistance are derived from PI88788. Resistance derived from PI88788 is breaking down due to narrow genetic background and SCN population shift. PI88788 requires mainly the rhg1-b locus, while 'Peking' requires rhg1-a and Rhg4 for SCN resistance. In the present study, whole genome re-sequencing of 106 soybean lines was used to define the Rhg haplotypes and investigate their responses to the SCN HG-Types. The analysis showed a comprehensive profile of SNPs and copy number variations (CNV) at these loci. CNV of rhg1 (GmSNAP18) only contributed towards resistance in lines derived from PI88788 and 'Cloud'. At least 5.6 copies of the PI88788-type rhg1 were required to confer SCN resistance, regardless of the Rhg4 (GmSHMT08) haplotype. However, when the GmSNAP18 copies dropped below 5.6, a 'Peking'-type GmSHMT08 haplotype was required to ensure SCN resistance. This points to a novel mechanism of epistasis between GmSNAP18 and GmSHMT08 involving minimum requirements for copy number. The presence of more Rhg4 copies confers resistance to multiple SCN races. Moreover, transcript abundance of the GmSHMT08 in root tissue correlates with more copies of the Rhg4 locus, reinforcing SCN resistance. Finally, haplotype analysis of the GmSHMT08 and GmSNAP18 promoters inferred additional levels of the resistance mechanism. This is the first report revealing the genetic basis of broad-based resistance to SCN and providing new insight into epistasis, haplotype-compatibility, CNV, promoter variation and its impact on broad-based disease resistance in plants.

Molecular characterization of genomic regions for resistance to Pythium ultimum var. ultimum in the soybean cultivar Magellan.

KEY MESSAGE: Two novel QTL for resistance to Pythium ultimum var. ultimum were identified in soybean using an Illumina SNP Chip and whole genome re-sequencing. Pythium ultimum var. ultimum is one of numerous Pythium spp. that causes severe pre- and post-emergence damping-off of seedlings and root rot of soybean [Glycine max (L.) Merr.]. The objective of this research was to identify quantitative trait loci (QTL) for resistance to P. ultimum var. ultimum in a recombinant inbred line population derived from a cross of 'Magellan' (moderately resistant) and PI 438489B (susceptible). Two different mapping approaches were utilized: the universal soybean linkage panel (USLP 1.0) and the bin map constructed from whole genome re-sequencing (WGRS) technology. Two genomic regions associated with variation in three disease-related parameters were detected using both approaches, with the bin map providing higher resolution. Using WGRS, the first QTL were mapped within a 350-kbp region on Chr. 6 and explained 7.5-13.5% of the phenotypic variance. The second QTL were positioned in a 260-kbp confidence interval on Chr. 8 and explained 6.3-16.8% of the phenotypic variation. Candidate genes potentially associated with disease resistance were proposed. High-resolution genetic linkage maps with a number of significant SNP markers could benefit marker-assisted breeding and dissection of the molecular mechanisms underlying soybean resistance to Pythium damping-off in 'Magellan.' Additionally, the outputs of this study may encourage more screening of diverse soybean germplasm and utilization of genome-wide association studies to understand the genetic basis of quantitative disease resistance.